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Preface 


The  purpose  of  this  study  was  to  develop  the  detection  characteristics  of  a 
steady  target  in  Weibull  radar  clutter  and  to  provide  a  means  of  predicting  the  radar 
detection  range  of  such  targets.  This  problem  originated  after  first  investigating  the 
problem  of  predicting  the  radar  detection  range  of  helicopters  flying  nap  of  the  earth 
profiles.  This  study  can  be  considered  as  a  first  step  in  the  approach  to  such  a 
problem. 

I  became  interested  in  radar  detection  problems  through  my  courses  at  AFIT 
and  chose  this  particular  research  project  after  being  introduced  to  it  by  my  thesis 
advisor. 

I  wish  to  express  my  gratitude  to  everyone  who  helped  me  to  complete  this 
thesis.  I  thank  my  thesis  advisor,  Dr.  Vittal  Pyati,  for  introducing  me  to  this 
topic  and  for  the  guidance  he  provided  through  the  course  of  my  research.  I  am 
also  grateful  to  the  members  of  my  thesis  committee,  Lt.  Col.  David  Meer,  Lt.  Col. 
Dave  Norman,  and  Dr.  B.N.  Nagarsenker,  for  their  comments  and  suggestions.  1  also 
wish  to  thank  my  classmates  in  GE-90D  and  especially  those  in  the  Communications^ 
Radar'  section  for  tiieir  help  with  this  thesis  effort  as  well  as  our  course  work.  Finally, 
I  want  to  extend  a  special  thanks  to  my  wife,  Kathy,  for  the  encouragement,  support, 
and  love  that  she  provided  throughout  this  project. 

Leon  C.  Rountree 
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Abstract 

In  this  thesis,  the  detection  characteristics  of  a  steady  target  immersed  in 
radar  clutter  characterized  by  returns  statistically  modeled  according  to  a  Weibull 
probability  density  function  are  developed.  A  detector  model  is  developed  which  is 
a  modification  of  an  envelope  detection  scheme  used  to  detect  deterministic  signals 
in  Gaussian  noise.  By  examining  the  statistics  of  the  output  of  the  detector  under 
alternate  hypotheses,  an  expression  for  the  probability  of  detection  (Pd)  is  developed 
in  terms  of  the  probability  of  false  alarm  (P/o)>  the  signal-to-clutter  ratio  (SCR), 
and  the  Weibull  skewness  parameter  (b),  which  is  a  function  of  the  type  of  terrain 
causing  the  clutter.  This  expression  is  used  in  conjunction  with  a  form  of  the  radar 
range  equation  for  area  clutter  to  produce  equations  which  describe  the  relationship 
between  Pd,  Pjat  SCR,  b,  and  the  radar  detection  range-  These  equations  involve 
Marcum’s  Q  function  and  an  algorithm  is  presented  which  evaluates  this  function. 
The  algorithm  is  a  series  approximation  and  can  be  computed  to  a  desired  accuracy. 
The  algorithm  is  also  presented  as  a  recursive  relationship  and  it  is  implemented 
using  a  FORTRAN  computer  program.  Plots  of  Pd  versus  SCR  and  Pd  versus  range 
aie  developed  for  various  values  of  Pja  and  b.  The  plots  provide  a  moans  of  evaluating 
the  detector  performance  and  predicting  the  detection  range  far  several  values  of  b; 
thus  for  several  types  of  terrain. 


THE  RADAR  RANGE  EQUATION  FOR  THE  DETECTION 
OF  STEADY  TARGETS  IN  WEIBULL  CLUTTER 


L  Introduction 

The  practical  problem  of  predicting  the  radar  detection  range  of  targets  in 
radar  clutter  requires  a  modification  to  the  standard  radar  range  equation  found  in 
textbooks  and  other  literature.  The  standard  radar  range  equation  is  derived  under 
the  assumption  that  the  primary  factor  limiting  detectability  is  noise  created  within 
the  receivei'  and  captured  by  the  antenna.  This  assumption  is  inadequate  for  targets 
in  radar  clutter  because  detectability  is  not  limited  by  noise  alone.  It  is  limited  by 
cluttter  returns  caused  by  trees  and  terrain,  whicli  in  most  cases  overwhelms  receiver 
noise. 

Receiver  and  ajitcnna  noise  voltage  amplitude  is  typically  modeled  as  a  Gaus¬ 
sian  random  process  and  tlie  problem  of  the  detection  of  a  steady  target,  or  a  target 
diaiactcrized  by  a  deterministic  signal,  in  Gaussian  noise  is  well-documented.  How¬ 
ever,  the  detection  problem  of  targets  in  non-Gaussiau  interference  is  the  topic  of 
recent  researdi.  To  accurately  predict  the  radar'  detectiorr  range  of  steady  tiurgets  itr 
radar  clutter,  a  statistical  irrodel  of  the  clutter  must  be  developed.  Receirt  literature 
suggests  that  the  Weibull  probability  density  function  (pdf)  can  be  used  to  model 
the  rrature  of  radar  returns  from  laird  and  sea  clutter  (1:1)  (6:221)  (3:136). 

i.i  Airplicalion 

The  prediction  of  the  radar  detection  range  of  targets  hr  radar'  clutter  can  be 
applied  to  the  problem  of  dcterniiiiiiig  the  radar  detection  rmige  of  helicopters  (lying 
trap  of  the  earth  (NOE)  prolilcs.  Helicopters  pilots  use  NOE  flight  to  avoid  radar 
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detection.  The  helicopter  target  signal  becomes  immersed  in  the  clutter  returns 
caused  by  the  terrain  when  the  helicopter  flies  close  to  the  ground.  The  problem  of 
detecting  steady  targets  in  Weibull  clutter  can  be  used  as  a  first  step  in  the  analysis 
of  the  problem  of  detection  of  helicopters  flying  NOE  profiles. 

1.2  Problem  Statement 

This  thesis  addresses  the  problem  of  predicting  the  radar  detection  range  of 
steady  targets  immersed  in  radar  clutter  characterized  by  returns  that  can  be  sta¬ 
tistically  modeled  according  to  a  Weibull  pdf. 

1.3  Scope 

The  detection  characteristics  of  a  deterministic,  or  exactly  known,  signal  cor¬ 
rupted  by  clutter  interference  modeled  according  to  a  Weibull  pdf  will  be  developed. 
The  objectives  of  this  thesis  are  to  determine  the  relationship  between  the  probability 
of  detection,  the  probability  of  false  alarm,  and  the  signal-to-clutter  ratio  (SCR)  for 
deterministic  signals  detected  in  Weibull  clutter  and  to  develop  a  set  of  curves  that 
shows  the  minimum  required  SCR  and  the  maximum  detection  range  for  a  desired 
probability  of  detection  for  various  values  of  probability  of  false  alarm. 

The  curves  will  be  developed  for  diilerent  values  of  the  skewness  pai'araetcr, 
b,  of  the  Weibull  pdf  (see  Appendi.x  B).  Different  values  of  b  have  been  found  to 
corresimnd  to  different  types  of  terrain  when  the  radar*  returns  from  the  terrain  are 
statistically  modeled  according  to  a  Weibull  pdf  (18:737-738).  The  values  of  b  used 
to  develop  these  curves  will  spair  the  range  of  typical  values  representing  several 
tyjres  of  terrain. 

j.4  Ai^mnptiom 

The  prinicuy  ussuin|>tion  made  is  that  the  target  to  be  detected  is  a  steady 
target.  That  is,  the  antplilude  of  the  return  signal  from  the  target  is  exactly  kttown. 


Also,  the  detection  of  the  target  will  be  carried  out  with  a  single  pulse  observation 
of  the  target’s  return  signal. 

The  radar  clutter  will  be  assumed  to  be  area  clutter  spread  over  a  patch  of 
terrain.  It  will  also  be  assumed  that  the  interference  caused  by  the  rr  iar  clutter 
will  overwhelm  the  receiver  and  antenna  noise  so  that  the  noise  can  be  considered 
negligible  leaving  the  clutter  as  the  only  interference  corrupting  the  signal. 

1.5  General  Approach 

The  first  step  in  the  approach  to  this  problem  is  to  develop  a  detector  model 
that  simulates  a  detection  scliemc  for  steady  targets  in  Weibull  clutter.  This  model 
can  be  developed  by  modifying  a  detection  scheme  commonly  used  to  detect  steady 
tai'gets  in  Gaussian  noise. 

By  examining  the  statistics  of  the  output  of  the  detector  model  under  the  hy¬ 
potheses  of  a  signal  present  and  clutter  alone  present,  expressions  for  the  probability 
of  detection,  Pd,  and  the  probability  of  false  alarm.  P/a,  can  be  found  for  the  detec¬ 
tion  of  steady  targets  in  Weibull  clutter.  These  expressions  describe  the  relationship 
between  the  Pd,  the  P/a,  and  the  signal4o-clutter  ratio,  SCR. 

Titis  relationship  can  be  applied  to  a  modiried  form  of  the  radar  riurge  equation 
for  urea  clutter  to  determine  the  relationship  between  the  maximum  radar  detection 
range,  Pd,  and  P/a.  Finally,  the  urathematical  expressions  that  describe  these  rc- 
lationslii|>s  are  examined  using  digital  computation  in  order  to  develop  plots  of  Pd 
versus  SCR  and  P^  versus  range  for  different  values  of  P/a  and  Weibull  skewness 
l>arameters. 

1.0  Sequmoe  of  h^^eulation 

First  a  summary  of  how  the  Weibull  pdf  can  be  used  to  model  radar  surface 
clutter  is  presented  in  Chapter  2.  Chapter  2  also  includes  a  presentation  of  analyses 
of  the  detectiou  problem  of  steady  targets  in  Weibull  iulerterencc. 


The  expressions  relating  Pd,  P/a,  SCR,  and  range  are  developed  in  Chapter  3. 
Chapter  4  contains  the  evaluation  of  these  expression  using  a  FORTRAN  computer 
program.  Plots  displaying  the  relationship  between  Pd,  P/a,  SCR,  and  range  are  also 
included  in  Chapter  4. 

Finally,  conclusions  are  presented  and  recommendv.,«:ons  for  further  research  are 
suggested  in  Chapter  5.  An  introduction  to  the  problem  of  detecting  a  Weibull  target 
in  Weibull  clutter  is  included  in  the  recommendations.  A  listing  of  the  computer 
program  used  to  develop  the  plots  described  above  is  included  as  an  appendix. 

1.7  Equipment 

A  VAX  11/785  computer  with  a  FORTRAN  77  compiler  is  used  to  implement 
the  computer  program  found  in  Appendix  A. 


II.  Background 


In  this  chapter,  the  Weibull  clutter  model  is  introduced  and  some  existing 
analyses  of  steady  target  detection  in  Weibull  clutter  is  reviewed.  First,  the  detection 
problem  of  a  deterministic  signal  in  Gaussian  noise  is  recalled. 

2.1  Target  Detection  in  Gaussian  Noise 

Marcum  and  Swerling  produced  a  definitive  study  of  the  statistical  problem 
of  target  detection  by  a  pulse  radar.  Marcum’s  paper  dealt  with  steady  target 
detection  in  Gaussian  noise,  while  Swerling  extended  Marcum’s  results  to  the  case 
of  fluctuating  targets  in  Gaussian  noise,  Swerling  created  four  statistical  models 
for  different  types  of  fluctuating  tai'gets  which  ai'e  commonly  referred  to  in  radar 
literature  as  Cases  1,  2,  3,  and  4.  Marcum’s  work  has  become  known  as  Case 
0  (12:62-83,273-288).  Marcum’s  Case  0  will  be  the  problem  examined  in  this  thesis 
with  the  exception  of  the  signal  interference  no  longer  being  Gaussian  noise  but 
Weibull  clutter. 

Marcum  created  expressions  for  Pj  and  P/a  for  Case  0.  He  also  determined 
the  relationship  between  Pja^  and  the  maximum  detection  range  and  plotted  the 
results  (12:85-138).  Marcum’s  expression  for  Pd  has  become  known  as  Mai'cum’s  Q 
function  (14:40).  The  Q  function  will  be  discussed  more  in  Chapters  3  and  4. 


2.2  Weibull  Clutter  Model 


2.2,1  Justijkation  for  the  Weibull  Model.  Several  pdfs  have  been  used  to 
model  the  nature  of  radar  sigiud  returns  from  land  iutd  sea  clutter.  Both  log-normal 
atul  llayleigli  pdfs  have  been  uschI  to  model  clutter  statistically  in  terms  of  spatial 
distribution.  The  liaylcigh  pdf  ha.s  been  used  for  clutter  observed  at  large  gracing 
i>ngles  (greater  tlian  5  degrees).  The  log-normal  distribution  has  been  used  more  for 


small  grazing  angles  because  at  small  angles,  land  clutter  density  generally  varies 
over  a  wider  range  than  the  Rayleigh  pdf  is  suited  for  (1:1). 

In  Schleher’s  paper,  he  states  that  the  Rayleigh  model  usually  underestimates 
the  dynamic  range  of  the  clutter  distribution,  while  the  log-normal  model  overes¬ 
timates  it.  The  Weibull  pdf  can  be  made  to  approach  either  the  Rayleigh  or  the 
log-normal  by  an  adjustment  of  its  parameters  and  it  can  be  used  as  a  model  over  a 
wider  range  of  conditions.  Schleher  further  states: 

From  a  detection  standpoint,  it  can  be  said  that  the  log-normal  dis¬ 
tribution  represents  the  worst  case  distribution,  the  Rayleigh  the  best 
case,  and  the  Weibull  an  intermediate  model  that  may  more  accurately 
represent  the  real  detection  performance  in  clutter.  (18:736-737) 

The  Weibull  pdf  has  been  found  to  correspond  to  some  land  clutter  distribu¬ 
tions  quite  well.  There  are  several  sources  of  measured  distributions  of  clutter  that 
exhibit  a  Weibull  distribution.  Data  have  been  recorded  for  several  types  of  terrain 
including  grasslands,  forests,  cultivated  land,  and  mountains.  All  these  cases  were 
found  to  fit  the  Weibull  pdf  (1:1-12).  Data  taken  from  radars  illuminating  the  ocean 
indicate  that  sea  clutter  is  also  distributed  according  to  a  Weibull  pdf  (3:929). 

2.2.2  Weibull  Parameters  for  Clutter  Distributions.  A  graphical  technique 
for  determining  if  measured  distributions  of  clutter  returns  can  be  represented  by  a 
Weibull  pdf  is  available.  The  equation  describing  the  Weibull  distribution  can  be 
manipulated  so  that  the  variate  of  the  distribution  is  expressed  as  a  linear  function  of 
the  median  of  the  distribution  and  the  Weibull  skewness  parameter,  b.  The  slope  of 
the  line  created  by  this  function  is  equal  to  1/b.  Merrsured  clutter  distributions  can 
be  plotted  to  determine  if  the  distribution  fits  the  Weibull  model  and  to  determine 
the  value  of  b.  If  the  points  lie  in  a  reasonably  straight  line,  the  Weibull  distribution 
can  be  used  to  model  the  clutter  (1:2-4). 
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Table  1.  Weibull  Skewness  Parameter  for  Clutter  Returns  (15:234) 


Terrain/Sea  state 

Preq 

Beam- 

width(®) 

Grazing 

angle(®) 

Pulse-  b 

length(/rs) 

Rocky  mountains 

S 

1.5 

2 

.512 

Wooded  hills 

L 

1.7 

.5 

3 

.626 

Forest 

X 

1.4 

.7 

.17 

.506-.513 

Cultivated  land 

X 

1.4 

.7-5 

.17 

.606-2 

Sea  state  1 

X 

.5 

4.7 

.02 

1.452 

Sea  state  2 

Ku 

5 

1-30 

.1 

1.16-1.783 

Several  measured  distributions  of  clutter  returns  from  various  terrain  types 
which  exhibit  the  Weibull  distribution  have  been  obtained.  Table  1  shows  the  values 
of  b  which  correspond  to  different  types  of  land  and  sea  clutter  (15:233-234). 


2.3  Detection  in  Weibull  Clutter 

An  approximate  analysis  of  the  detection  characteristics  for  a  steady  target 
in  Weibull  clutter  was  made  by  Ekstrom.  In  this  analysis,  it  was  assumed  that  the 
signal  voltage  was  much  greater  than  the  median  voltage  value  of  the  clutter  return 
distribution.  Under  this  assumption,  closed  form  expressions  that  approximate  the 
Pd  when  Pd  is  greater  than  1/2  were  developed  (6:222-225). 

Schleher  has  examined  the  problem  more  thoroughly  by  using  the  Weibull  and 
the  Weibull-Riccan  pdfs  to  describe  the  statistics  of  the  output  of  an  envelope  de¬ 
tector  receiver  under  the  hypotheses  of  clutter  alone  present  and  signal  plus  clutter 
present  respectively.  Although  the  details  of  the  analysis  are  not  included  in  the  his 
report,  the  Pd  was  determined  ’’using  a  chaiacteristic  function  approach  in  conjunc¬ 
tion  with  the  fiist-Fourier  transform.”  The  results  of  this  study  are  plots  of  Pd  versus 
signal-to-cluttcr  ratio  for  a  Pfa  of  10"®  and  b  equal  to  2, 1.2,  .8,  and  .6  (18:737-739). 
These  plots  can  be  compared  to  the  ones  developed  in  this  thesis. 

Other  resear-clicrs  have  taken  the  problem  a  step  further  by  allowing  the  target 
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to  be  a  fluctuating  one.  Goldstein  obtained  expressions  for  the  Pd  of  a  Rayleigh 
target  in  Weibull  clutter.  These  results  are  based  on  the  assumption  that  the  signal- 
to-clutter  ratio  is  large.  For  fluctuating  targets,  the  pdf  of  the  envelope  of  the  received 
signal  when  both  target  and  clutter  are  present  can  be  approximated  by  the  pdf  of 
the  target  signal’s  envelope  when  the  signal-to-clutter  ratio  is  large  (9:90-91).  Chen 
and  Morchin  also  used  this  fact  in  developing  an  expression  for  Pd  of  a  Rayleigh 
target  in  Weibull  clutter.  Their  results  were  plotted  for  b  equal  to  .5,  .6,  .7,  and 
.8  (3:929-932). 
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III.  Development  of  Detection  Characteristics 

The  problem  of  detecting  a  signal  known  exactly  corrupted  by  Weibull  clut¬ 
ter  interference  can  be  approached  by  first  examining  the  detection  problem  of  a 
deterministic  signal  immersed  in  Gaussian  noise.  When  white  Gaussian  noise  alone 
is  passed  through  a  narrowband  filter,  the  probability  density  function  (pdf)  of  the 
envelope  of  the  output  noise  voltage  is  Rayleigh.  The  envelope  of  Gaussian  noise  plus 
a  deterministic  signal  passed  through  a  narrowband  filter  is  distributed  according  to 
a  Ricean  pdf  (19:23-27).  These  results  will  be  useful  in  the  analysis  of  the  detection 
problem  in  Weibull  clutter. 

3,1  Rayleigh-to-WeibuU  Transformation 

A  Rayleigh  random  variable,  r,  representing  the  envelope  of  a  narrowband 
Gaussian  process  can  be  transformed  into  a  Weibull  random  variable,  w,  by  raising 
r  to  a  power,  v  (8:10-11): 

w  =1  r”  I  (1) 

The  value  of  v  that  will  result  in  w  being  a  Weibull  random  variable  is  found  by 
performing  the  transformation  of  equation  1  on  r.  The  variable,  r,  has  a  Rayleigh 
pdf  as  shown: 

p(r)  =  (?'/(7^)exp(-)*^/2(7^)  r>0  (2) 

The  pai’ameter,  is  the  mean-square  voltage  of  the  Gaussian  input.  The  pdf  of  w 
is  found  as  follows 

piw)  =  p{r  ss  tn*/*’)  I  dr/dw  \ 

=  (u,(2/.'-i)/^^3)exp(-ty2/72(r“)  ta>0  (3) 

A  Weibull  pdf  has  the  form 
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Figure  1,  Detector  for  Deterministic  Signal  in  Gaussian  Noise 


p(aj)  =  (6/af)a?^“^exp(-aj7o)  a;  >  0  (4) 


From  equations  3  and  4,  it  can  be  seen  that  w  is  Weibull  when 


j/  =  2/6 


(5) 


3,2  Detector  Model 

The  Weibull  envelope  clutter  process  can  be  modeled  artificially  by  modifying 
a  detection  scheme  used  to  detect  a  signal  known  exactly  with  unknown  phase  in 
Gaussian  noise,  A  simplified  block  diagram  of  a  detector  to  detect  exactly  known 
signals  with  unknown  phase  in  Gaussian  noise  is  shown  in  Figure  1  (5:298-302). 


signal 
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Figure  2.  Detector  Model  for  Deterministic  Signal  in  Weibull  Clutter 
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w(t) 

Threshold 

b 

Device 

Using  the  results  of  the  previous  section,  it  is  seen  that  a  Weibull  process 
can  be  simulated  by  introducing  an  operation  after  the  envelope  detector  that  will 
transform  the  the  Rayleigh  variable,  r,  into  a  Weibull  variable,  w,  when  noise  alone 
is  present.  The  operation  necessary  is  determined  from  equations  1  and  5  and  is 
given  by 

w  =  (6) 

A  block  diagram  for  this  detector  is  shown  in  Figure  2. 

Use  of  this  model  allows  one  to  approach  the  detection  problem  of  a  deter¬ 
ministic  signal  in  Weibull  clutter  by  first  taldng  advantage  of  the  well-documented 
results  of  the  detection  of  a  deterministic  signal  in  Gaussian  noise. 
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3.2  Radar  Parameter  Calculations 

3.3.1  Probability  of  False  Alarm.  When  there  is  no  target  signal  present, 
the  signals,  r(t)  and  w(t),  from  Figure  2  are  distributed  according  to  Rayleigh  and 
Weibull  pdfs  respectively.  The  pdf  of  w(t)  is 

p{wfHo)  =  {b{a)w’^~^exp{—w^la)  w>0  (7) 

where  Ho  represents  the  hypothesis  that  no  target  is  present.  The  parameter,  b, 
indicates  the  degree  of  skewness  in  the  distribution  and  a  is  related  to  the  median 
of  the  distribution  (15:234).  While  b  is  determined  by  the  type  of  terrain  causing 
the  clutter  returns,  (x  is  determined  by  the  mean-square  voltage,  of  the  Gaussian 
input.  It  can  be  seen  from  equations  3  and  4  that  a  =  2a^.  The  parameter,  a,  is 
related  to  the  median  by  (1:2) 

twm  =  (aln(2))^^^  (8) 

where  Wm  is  the  median.  From  equation  8  and  the  relation  that  a  =  2o-^,  it  is  seen 

that 

a»  =  »S./(21n(2))  (9) 

The  probability  of  false  alar  m  is  determined  by  integrating  the  pdf  of  equation  7 
from  the  detection  threshold,  wj',  to  infinity: 

i  Qxp{-‘W^fa)dw 

=  exp(-tt>5./2)  (10) 

When  P/a  is  a  fixed  parameter,  tuj’  can  be  solved  for: 

ti;r  =  (21n(l/P/a))V‘»  (11) 
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3.3.2  Probability  of  Detection.  When  an  exactly  known  signal  is  present  at 
the  input  of  the  detector,  r(t)  has  a  Ricean  pdf.  Using  DiPranco  and  Rubin’s  nota¬ 
tion,  the  Ricean  pdf  is 

p{rJHi)  =  r exp[-(r^  -f  1l)/2]Io{r\/iZ)  r  >  0  (12) 

where  Hi  is  the  hypothesis  that  a  signal  is  present  and  7o()  represents  the  zero-order 
modified  Bessel  function  of  the  first  kind.  %  is  defined  to  be  2E/iVo  where  E  is  the 
detected  target  signal  energy  and  No  is  the  one-sided  power  spectral  density  of  a  zero- 
mean  white  Gaussian  noise  input.  DiFranco  and  Rubin  refer  to  the  parameter,  as 
the  peak  signal-to-noise  power  ratio,  which  is  the  ratio  of  maximum  instantaneous 
signal  power  to  average  noise  power  out  of  a  matched  filter.  Tl  appears  often  in  radar 
literature  (5:302-307). 

When  a  signal  is  present,  r(t)  undergoes  the  transformation  of  equation  6  and 
w(t)  has  the  pdf  developed  below: 

p{wllli)  =  p{r  =  I  dr/dw  | 

=  exp[-(M>’®^*'  -I-  'E)/2]/o(u;‘^*'\/^)(l/t/)t/;^*^‘'^“' 

~  72)/2]/<,(u;*^‘'v^)  iw>0  (13) 

Substituting  v  =  2/b,  p{vtlH\)  becomes 

p{wJH\)  =  (6/2)u;^“^  oxp[--(w^  -I-  'R)l2]lo{'w'*^^y/K)  w>0  (14) 

To  determine  the  probability  of  detection  for  a  signal  known  exactly  in  Woibull 
envelope  clutter,  equation  14  is  integrated  from  the  detection  threshold,  u>2',  to 
infinity: 

P^ss  I  (6/2)u;*''*  exp(— (toN' 7il)/2]/«(iw^^®\/^)du)  (15) 
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This  equation  is  of  the  form  of  Marcum’s  Q  function  given  by  (14:40) 


/•OO 

Q{ctt^)=  vexp[-{v^ +  a^)/2]Io{av)dv  (16) 

Making  the  following  substitutions: 

a  = 

^  =  wt 

V  = 

dv  =  {bl2)w^^^^^~^dw 

the  integral  of  equation  15  can  be  written  as 

Pd  —  f  VGxp[-~{v^  •\-'Jl)/2]lo{vy/^)dv  (17) 

=  Q{y/n,WT)  (18) 

Combining  equations  11  and  18,  Pd  becomes 

/'d  =  «(\/K.(21n(l/J>/.))‘'‘)  (19) 

Also,  Marcum’s  Q  function  can  bo  written  in  terms  of  the  incomplete  Toronto  func¬ 
tion  as  (12:160) 

Q(a,;3)  =  l-2i;^(l,0,ayv^)  (20) 

Therefore,  Pd  becomes 

ft  =  l-W«■'‘/^/3(^.0.(K/2)■'’)  (21) 

When  b  =  2,  it  is  seen  that  the  pdf  of  w(t)  is  llayleigh  and  the  Pd  and  Pj^  are 
the  same  as  those  found  for  the  case  of  an  exactly  known  signal  with  unknown  phase 
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detected  in  the  presence  of  Gaussian  noise  with  an  envelope  detector.  Several  texts, 
including  DiPranco  and  Rubin’s  contain  plots  that  show  the  required  peak  signal- to- 
noise  ratio,  %,  needed  to  acheive  a  desired  probability  of  detection  for  various  values 
of  false  alarm  (5:308).  These  plots  are  identical  to  the  ones  that  would  be  created 
by  an  evaluation  of  equation  19  with  b  =  2. 


3.4  Relationship  between  SCR,  Pja,  and  Pd 

In  the  case  of  the  detection  of  a  deterministic  signal  corrupted  by  Weibull 
clutter,  b  is  typically  not  equal  to  2  and  Pd  is  determined  by  equation  19  or  21. 
When  the  target  signal  is  corrupted  by  Weibull  clutter,  it  is  more  meaningful  to  use 
a  signal-to-clutter  ratio  (SCR),  instead  of  the  parameter,  Ti. 

The  clutter  power  is  determined  from  the  second  moment  of  the  Weibull  pdf 
of  equation  7  as  (15:34) 


E[w^]  ==  r(l  +  (2/6)) 


(ln(2))2/6 


where  r()  is  the  Gamma  function  and  ta„,  is  the  median  of  tlie  Weibull  distribution. 
If  the  amplitude  of  the  signal  before  the  envelope  detector  is  A,  the  power  signal-to- 
clutter  ratio  at  point  b  of  Figure  2  is 


(yl3/6)2 

K/(ii.2)W‘)r(i  +  (2/t)) 
<r(i+(2/6)) 


Recall  that  is  the  peak  signal-to-noise  power  ratio.  The  |)eak  signal  |)ower 
is  related  to  the  average  instantaneous  iM>wer  by  (5:296) 


{l>eak  instantaneous  signal  2>ower)  2{average  instantaneous  signal  j/ower) 
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The  average  instantaneous  signal  power  of  a  sinusoidal  signal  of  amplitude,  A, 
is  i4^/2  and  the  average  noise  power  is  a^.  Therefore,  TZ  expressed  in  terms  of  is 

n  =  [2(AV2)]/a" 

=  (25) 


Using  equation  9,  it  is  seen  that 


7^=(A221n2)/i«i, 


(26) 


Therefore,  SCRi)  is  related  to  Tl  by 


SCRt>  = 


(7^/2)2/^ 

r(i+(2/6)) 


(27) 


To  determine  the  SCR  at  the  input  of  the  envelope  detector,  the  inverse  of  the 
transformation  of  equation  6  can  be  performed  on  SCRi,  to  give 


SCR  =  {SCRi,f^^ 

Ti 

”  2(r(H- (276))]^/2 

This  leads  to 


(28) 


71  =:  2iSGR)[r{l-{-{2lb))ff^  (29) 

\/^  =  (30) 

Combining  equations  18  and  30,  the  relationship  between  SCR,  and  Pd  is 

i'j  =  0(l2(SC«)l‘7=ll'(l  +  (2/6))l‘'M2lii(l/i'/.)l'7‘)  (31) 

where  Q(a,/8)  is  Maicuiu’s  Q  (unction. 
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3.5  Radar  Range  Equation 


The  radar  range  equation  for  targets  in  area  clutter  can  be  divided  into  two 
cases:  beamwidth-limited  case  and  pulse-length-limited  case.  For  large  grazing  an¬ 
gles,  the  area  of  terrain  illuminated  by  the  radar  is  limited  by  antenna  beamwidth; 
while  for  small  grazing  angles,  the  area  is  limited  by  the  pulse-length  of  the  trans¬ 
mitted  signal.  The  radar  range  equations  for  each  of  these  two  cases  are  (16:63-67) 


and 


where 


i(sin 

{irj4)d<fi{SCR)ao 

L{cOS  lf))<Tt 
iSGR){cTl2)6ao 


for  tan  0  > 


for  tan 


(f>R 

cr/2 


(l>R 

ct/2 


(32) 


(33) 


•  I.  =5  all  of  the  various  losses  encountered  by  the  radar 


9-  ilf  =  grazing  angle  of  the  radar 

•  »  target  radar  cross  section 

•  SCR  sa  mean  sigoal'torclutter  power  ratio 

•  <7o  s=s  the  backscatter  (clutter)  radar  aoss  section  per  unit  are^of  reflecting 

•  sa  3'dB  beamwidth  in  azimuth 
s  ^HilE  beantwidth  ill  elevation 


•  r  8  pdlso  duration  of  the  transmitted  signal 
•:  e  8  spijt'd  of  light 
t  R  S3  raclar>  detection  range 

The>X1l.<^i  be  solved  for  to  obtain 


imip  > 


m 


(34) 


and 


SCR  = 


L{cos 


for  tan  ^  < 


(f>R 


(35) 


R(arl2)6(T,  '  crft 

Now,  define  a  normalizing  range,  Ro,  which  is  the  range  at  which  the  SCR  becomes 
unity  (2:23-27).  Using  equations  34  and  35,  Ro  becomes 


tanV>  > 


(f>R 

cr/2 


(36) 


and 


i(cos^)crt  <j)R 

M  tanV>< 


(cr/2)(lff.  •“’''ct/2 

Now,  all  radar  parameters  are  included  in  R^  and  SCR  can  be  related  to  range  readily 
by 


and 


RlRo  =  {llSCRff^  for  tan^> 


R/Ro  =  1/SCR  for  tan0  < 


<f>R 

ct/2 

<I>R 


ct/2 


(38) 


(39) 


Therefore,  for  small  grazing  angles  (pulse-length-limited  case),  the  detection 
range  is  inversely  proportional  to  the  SCR,  For  larger  grazing  angles  (beamwidth- 
limited  case),  the  detection  range  is  inversely  proportional  to  the  square  root  of  the 
SCR. 

The  relationship  between  detection  range,  J%  and  JP/4  can  be  determined  by 
combining  equations  31,  38,  and  39; 

where 


and 


P,  =  Q{{2Ro/RY^^[T{1  +  {2/b))ff\ [21n(l/Py«)]^/^)  for  tmij;  <  ^  (42) 


where 


Ro  = 


L{cos  il))at 
{cTf2)9ao 


(43) 


Equations  40  and  42  along  with  equation  31  can  be  evaluated  to  obtain  plots 
of  Pd  versus  R/P<,  and  Pd  versus  SCR  for  various  values  of  Pja  and  b.  The  detection 
problem  becomes  a  problem  of  evaluating  the  Q  function  for  the  arguments  shown 
in  equations  31,  40,  and  42.  This  is  the  subject  of  Chapter  4. 


IV.  Computer  Implementation 


The  expressions  for  Pi  from  equations  31,  40,  and  42  describe  the  detection 
characteristics  of  a  steady  target  in  Weibull  clutter.  In  order  to  evaluate  these 
expressions,  a  means  of  tabulating  Marcum’s  Q  function  must  be  developed. 

4.1  Algorithm  to  Compute  Q  function 

In  a  Johns  Hopkins  University  report  by  Fehlner,  an  algorithm  for  computing 
Marcum’s  Q  function  was  developed.  The  integrand  of  the  Q  function  was  red¬ 
erived  using  the  characteristic  function  and  the  Q  function  was  finally  expressed  as 
a  series  (7:23). 

The  characteristic  function  for  Marcum’s  Case  0  is  the  Fourier  transform  of 
the  integrand  of  the  Q  function  of  equation  16.  First,  the  substitution,  u  =  v^/2  and 
du  =  vdu  is  made  to  obtain  the  integrand,  I: 


I  =  exp[-(tt  +  l2)]lo{ay/^)  u  >  0  (44) 

The  diai'acteristic  function  becomes 

C  =  jf'”exp[-(u-haV2)]J,(aV^)e^‘^du  (45) 

Letting  x  ss  a^/2  and  p  =  ju;,  C  becomes 

Cs=  /  exp[--(u  4- a;)]/o(2>/Su)e’”Vu  (46) 

JQ 

Marcum  found  this  integral  to  be  (12:163-164) 

C='.(l/(p-{*l)]e"‘*oxp(a;/(p -1-1)1  (47) 
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Fehlner  takes  the  inverse  Fourier  transform  by  contour  integration  and  obtains 


f{t)  =  {l/2irj)  r 

%/— JOO 


e  ®  exp[a;/ (p  + 1)]  exp[(p  +  1  -  l)t] 


p  +  1 


dp 


e~^e  *  exp[-x/ (p  +  1)]  exp[t(p  +  1)]  , 
27ri  7-ioo  p  +  1  ^ 


(48) 


Fehlner’s  evaluation  of  this  integral  results  in 


*=o 

Now,  to  determine  the  Q  function,  the  integral 


(50) 


must  be  evaluated,  where  ut  is  the  threshold  corresponding  to  from  the  original 
expression  for  the  Q  function.  Upon  evaluation  of  this  integral,  Fehlner’s  final  result 
becomes  a  series  expression  given  by  (7:24-25) 

Q(a;,ur)  =  (51) 

i=0  j=0 


Recalling  that  x  =  a^/2,  u  =  v^/2,  and  thus  ut  =  /9*/2,  equation  51  becomes 

••  j* 


(52) 


»=0j=:0 


Tins  series  can  be  cairicd  out  to  the  desired  precision. 


4.2  Computational  Problems 

It  can  be  seen  from  equation  19  that  as  b  and  Pja  decreased,  the  second 
argument  of  the  Q  function,  fi,  becomes  large.  Also,  as  the  SCR  is  increased,  a,  the 
first  argument  of  the  Q  function,  becomes  largo.  Large  values  of  a  and  ^  become 

"  '  21  . 


computationally  hard  to  handle  in  equation  52.  The  expression 

Y 


(53) 


must  be  evaluated  with  a  =  or  a  =  ;0^/2.  Large  values  of  a,  and  k  cause 
numbers  for  X  well  outside  the  limitations  of  a  computer. 

This  problem  can  be  handled  by  first  evaluating 

h 

InX  =  — a  +  ^Ina  —  ^Ini  (54) 

i=l 

This  expression  can  be  evaluated  for  large  values  of  a  and  k  (7:30-31).  Then 


X  =  exp(ln  X)  (55) 


Equation  52  becomes 

00  i 

tsO  JbO 

where 

In  Ai  =  —a^/2  +  i  ln(o(^/2)  -  ^  In  & 

k=l 

InBj  =  -/?V2  +  iln(i3V2)“5^1n& 

k-l 

Ai  =  cxp(ln  /i,) 

Bj  =  exp(ln  Bj) 


(56) 


(57) 

(58) 

(59) 

(60) 


Equation  56  can  also  be  expressed  as  a  recursive  relationship  where 


(61) 
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The  initial  value  of  q  is  computed  as  follows 

In  Ao  =  -a^/2  (62) 

InB,  =  -^^2  (63) 

Ao  =  exp(lnAo)  (64) 

Bo  =  exp(lnjBo)  (65) 

qo  =  AoBo  (66) 

The  other  values  of  q  are  computed  from 

InAi  =  lnj4,_i+ln(a^/2)-lnz  (67) 

InBi  =  ln5,_i  +  ln(|3^/2)  —  Ini  (68) 

At  =  exp(lni4i)  (69) 

Bi  =  exp(ln5<)  (70) 

«<  =  qi-i  +  AiBi  (71) 


The  value  of  i  in  equation  61  must  be  large  enough  to  cause  Q  to  converge  to 
the  desired  accuracy.  To  obtain  a  six  digit  acurracy,  i  must  be  large  enough  that 

>  .999999  (72) 

4.3  Computer  Program  Design 

The  FORTRAN  computer  program  found  in  Appendix  A  was  developed  to 
produce  data  that  can  be  used  to  plot  Pj  versus  SCR  and  P<<  versus  R/P® 
different  values  of  b  and  P/o.  The  input  to  the  program  includes:  the  value  of  b, 
the  value  of  P/a,  and  whether  the  grazing  angle  is  iai'ge(beamwidthdimited  case)  or 
smail(pulse-length-limitcd  case). 

The  program  evaluates  equation  19,  whicli  is  Pd  expressed  as  the  Q  function 
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of  (72.)^/^  and  [2ln(l/P/a)]^^^*  First  a  and  ^  are  computed  according  to  equation  19. 
Then  equations  62  through  66  are  used  to  determine  qo.  Then,  a  DO  loop  is  carried 
put  to  determine  qi  according  to  equations  67  through  71.  The  loop  is  exited  when 
the  inequality  72  is  satisfied.  Next,  1Z  is  converted  to  SCR  according  to  equation  29. 
This  conversion  involves  the  Gamma  function  and  the  GAMMA  function  provided 
by  the  International  Mathematical  and  Statistical  Libraries  (IMSL)  is  utilized  in  the 
program  (10:GAMMA-1).  R/Ro  is  then  computed  according  to  equation  38  or  39. 
The  P(i,  SCR,  and  R/Ro  are  output  when  Pd  is  between  .0001  and  .9999.  Finally, 
Tl  is  incremented  and  the  computations  are  repeated  until  Pd  becomes  greater  than 
.9999.  The  flow  chart  describing  the  program  is  shown  in  Figure  3. 

Figures  4  through  12  are  plots  that  show  the  required  SCR  needed  to  produce 
a  desired  Pd  for  various  values  of  b  and  Pfa.  Figures  13  through  21  are  plots  of  Pd 
versus  R/Rq  for  large  grazing  angles  and  Figures  22  through  30  are  the  same  plots 
for  small  grazing  angles. 

As  P/a  and  b  become  small,  the  number  of  iterations  required  to  cause  the 
summation  of  equation  56  to  converge  becomes  extremely  large.  The  computation 
time  required  to  evaluate  the  algorithm  for  small  values  of  b  and  P/a  also  becomes 
extremely  large.  This  is  why  in  Figures  11, 12,  20,  21,  29,  aiid  30,  the  number  of  val¬ 
ues  of  P/a  ai'e  limited.  The  skewness  parameter,  b,  is  reduced  to  .5  in  Figures  12, 21, 
and  30.  Values  of  b  smaller  than  .506  are  not  documented. 
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Figure  3.  Program  Flowcliai't 
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SCR((1B) 

Figure  5.  Pd  versus  SCR  for  b  =  1.8, 

P/a  =  10-*,  10“S  10-®,  10-8,  IQ-IO  ^ 

10“^*,  IQ-'^  (left  to  right) 


27 


-2  0  2  4  0  8  10  12  14  10  18  20 

SCR(tlB) 


Figure  6.  Pd  versus  SCll  for  b  ~  1.6, 

P/a  ==  10"^  10-^  10“®,  10"8,  10“*®, 
10"*^,  10“***  (left  to  right) 
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Pd  versus  SCll  tor  b 
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Figure  8.  versus  SCR  for  b  =  1.2, 

^  10"^  l0-\  l0“^  10”»,  10"« 

10“*^  10"'**  (left  to  right) 
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SCll(dD),' 

Figure  U.  Pd  versus  SCil  for  b  =  .6, 

/%  =  10“^  10-\  10“^  10~^  10“^^ 

(left  to  right) 
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Figure  15.  Pd  versus  R/Ro  for  b  =  1.6  and  large  grazing  angles, 

Pu  =  10-^  10"^  10-«,  10-»,  10-»“,  10-^^  10”»^  (right  to  1 
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.9999 


Figure  16.  Pd  versus  R/Ro  for  b  =  1.4  and  large  grazing  angles, 

^  10-2, 10-4^  10-6^  10-»,  10-1“,  10-12^  10-14  (rigijt  to  left) 


Figure  18.  Pd  versus  R/Ro  for  b  =  1  and  large  grazing  angles, 

Pja  =  ^0-^  10-S  10-®,  10-8,  iQ-io^  io-«,  10-»4  (right  to  le 
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Figure  19.  Pj  versus  R/Ro  for  b  =  .8 
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Figure  24.  Pd  versus  11/ A,  for  b  js  1.6  and  small  graaiiig  migltjs, 


P/a  ~  i0“2, 10-“,  10-«,  10-»,  l0“»^  10“«  (right  to  left) 


^’iguie  25,  Pd  vcrsuf5  R//4  for  b  s  1.4  and  small  graaiiig  angles, 

Pja^  10“3,  lr^  10“«,  10”»,  10“*S  10-w  (right  to  left) 
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Figure  27.  Pd  versus  U/jRo  for  b  =  1  iuid  small  grazing  angles, 

Pj<x  =  10-’,  10-^  10-®,  10-®,  10-*®,  10-*’,  10-*^  (right  to  left) 


4-4  Results 

44 -  I  Validation,  In  order  to  determine  if  the  curves  presented  axe  valid,  the 
curves  developed  for  b  =  2  can  be  compared  to  those  presented  by  other  authors 
who  analyzed  the  detection  problem  of  a  steady  target  in  Gaussian  noise.  When  b 
=  2,  the  pdf  of  the  envelope  of  the  interference  is  Rayleigh  as  seen  from  equation 
7,  and  Rayleigh  envelope  interference  corresponds  to  Gaussian  noise  at  the  input 
to  an  envelope  detector.  Therefore,  the  detection  problem  for  a  steady  target  in 
Weibull  clutter  with  b  =  2  is  identical  to  the  problem  of  detecting  a  steady  target 
in  Gaussian  noise.  DiPranco  and  Rubin’s  curves  of  Pd  versus  the  peak  signal-to- 
noise  ratio  for  the  Gaussian  detection  problem  are  identical  to  the  curves  shown  in 
Figure  4  except  for  a  3  dB  shift  in  the  peak  signal- to-noise  ratio  axis  (5:308).  This 
discrepancy  is  explained  by  noting  that  the  curves  in  Figure  4  are  plotted  with  the 
mean  signal-to-clutter  ratio,  SCR,  on  the  x-axis  while  DiFranco  and  Rubin  use  the 
peak  signal-to-noise  ratio,  H.  Prom  equation  28,  it  can  be  seen  that  when  b  -  2, 
SCR  =  71/2,  thus  the  3  dB  discrepancy. 

To  further  validate  the  curves,  they  can  be  compared  to  the  results  obtained  by 
Schleher.  Schleher  developed  plots  of  Pd  versus  SCR  for  a  steady  target  in  Weibull 
clutter  for  P/o  =  10"®  and  b  =  2, 1.2,  .8,  and  .6  (18:737-739).  These  plots  correspond 
to  the  curves  shown  in  Figures  4,  8, 10,  and  11  with  P/o  =  10”®  and  the  curves  in 
these  hgures  match  the  ones  developed  by  Schleher. 

44’^  Compaiison  of  Results  for  Various  Values  of  b.  Fiom  the  results  ob¬ 
tained,  it  is  seen  that  as  b  is  decreased,  or  the  skewness  of  the  clutter  distribution 
is  decreased,  the  amount  of  SCR  required  to  produce  a  desired  Pj  is  increased  and 
thus  the  detection  range  is  decreased.  Table  2  shows  the  additional  SCR  required  to 
produce  a  Pj  =  .9  with  P/a  =  10"®  when  b  is  not  equal  to  2  comparerl  to  the  Rayleigh 
envelope  case  when  b  2.  The  decrease  in  range  for  large  and  small  grazing  angles  is 
also  shown.  Table  2  shows  that  as  the  skewness  of  the  clutter  distribution  is  dianged 

■53'  .. 


Table  2.  Additional  SCR  and  Decrease  in  Detection  Range 


b 

Additional 

SCR(dB) 

Decrease  in  R/Ro  for 
Small  Gra.zing  Angles(dB) 

Decrease  in  R/Ro  for 
Large  Grazing  Angles(dB) 

1.8 

1.2 

1.2 

.6 

1.6 

2.6 

2.6 

1.3 

1.4 

4.6 

4.6 

2.3 

1.4 

4.6 

4.6 

2.3 

1.2 

7.4 

7.4 

3.7 

1 

11.6 

11.6 

5.8 

.8 

17.8 

17.8 

8.9 

.6 

28.6 

28.6 

14.3 

.5 

38 

38 
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(b  decreases),  the  performance  of  the  envelope  detector  receiver  is  degraded. 

As  an  example  of  the  performance  degradation,  one  can  compare  the  detection 
characteristics  found  for  two  different  types  of  terrain.  Suppose  b  =  .626  for  terrain 
characterized  by  wooded  hills  and  b  =  1.452  for  the  ocean  (see  Table  1).  Using  the 
algorithm  and  computer  program  described  previously,  it  was  found  that  approxi* 
mately  24.1  dB  more  SCR  is  needed  to  detect  a  steady  target  over  wooded  hills  than 
that  needed  to  detect  a  steady  target  over  the  ocean.  For  small  grazing  angles,  the 
detection  range  is  also  24.1  dB  less  for  detection  over  wooded  hills  than  the  ocean 
and  for  large  grazing  angles,  the  detection  range  is  decreased  by  12.03  dB. 


Y.  Conclusions  and  Recommendations 


5.1  Conclusions 

The  purpose  of  this  thesis  was  to  develop  a  means  of  predicting  the  radar 
detection  range  of  steady  targets  in  Weibull  clutter.  This  was  accomplished  by  mod¬ 
ifying  a  detection  scheme  used  to  detect  steady  targets  in  Gaussian  noise  and  then 
determining  the  relationship  between  Pa,  P/o,  SCR,  and  b  and  finally  applying  this 
relationship  to  a  modified  form  of  the  radar  range  equation.  The  expressions  describ¬ 
ing  the  relationship  between  Pa,  P/o»  SCR,  b,  and  range  were  found  to  be  in  the  form 
of  Marcum’s  Q  function  (see  equations  31,  40,  and  42).  Computer  implementation 
of  the  recursive  algorithm  presented  in  Chapter  4  used  for  the  evaluation  of  the  Q 
function  was  capable  of  producing  data  used  to  plot  Pa  versus  SCR  and  Pa  versus 
range  for  various  values  of  P/a  and  b.  Six  digit  precision  was  used  in  developing  the 
plots  but  the  algorithm  is  capable  of  providing  greater  accuracy. 

The  plots  of  Pa  versus  SCR  developed  in  Chapter  4  describe  the  detection 
characteristics  of  a  steady  target  in  Weibull  clutter  and  show  that  the  performance 
of  the  detector  is  severely  degraded  as  the  skewness  of  the  clutter  distribution  is 
decreased.  That  is,  as  b  is  decreased,  the  SCR  required  to  produce  desired  values  of 
Pa  and  P/a  becomes  lai’ge  (see  Table  2). 

Finally,  the  plots  of  Pa  versus  R/Ro  provide  a  means  of  predicting  the  radai' 
detection  range  for  pai  ticulai'  types  of  clutter  (values  of  b)  and  desired  values  of  Pa 
and  P/a.  Several  improvemeiits  can  be  made  to  this  study  and  some  are  suggested 
in  the  next  section. 

5.2  Recommendations  for  l^tum  Jleseai'ch 

5.2.1  Weibull  Taryet.  The  next  step  in  the  Weibull  clutter  detection  problem 
is  to  allow  the  target  to  lluctuate.  The  Weibull  pdf  can  be  used  to  model  the  target 


just  as  it  was  used  to  model  the  radar  clutter.  Recall  that  the  Weibull  target  model 
is  equivalent  to  the  often  used  Rayleigh  model  when  b  =  2.  The  Weibull  target 
model  is  better  suited  than  the  steady  target  model  as  an  approach  to  problems  like 
detecting  helicopters  flying  nap  of  the  earth  profiles.  The  following  development  is 
provided  as  a  suggestion  on  how  to  approach  the  analysis  of  Weibull  target  detection 
in  Weibull  clutter. 

With  the  problem  of  detecting  a  Weibull  target  in  Weibull  clutter,  the  pdf  of 
the  detected  signal  under  the  hypothesis  of  clutter  alone  present  is  still  given  by 
equation  7.  When  a  target  is  immersed  in  clutter,  the  pdf  of  the  detected  signal 
describes  the  distribution  of  the  sum  of  two  random  processes.  The  pdf  of  the  sum 
of  two  random  variables  is  the  convolution  of  the  individual  pdfs  (4:186*189).  The 


pdf  corresponding  to  the  clutter  is 

Pc(w)  =  (6o/ao)tu^®“^  exp(--tt>*’*/ao)  u>  >  0  (73) 

and  the  pdf  corresponding  to  the  target  is 

Pt(w)  =  /at)  w>0  (74) 

so  the  pdf  of  the  signal  under  the  hypothesis  that  a  target  is  present  becomes 

p(w/ffi)==pc(w)*pt(w)  (75)  . 

Therefore  the  pdfs  of  the  detected  signal  under  the  alternate  hypotheses  are 

I. 

p{u)/Jio)  »  (5<,/ao)w;^“''‘ exp(~tn^*/®o)  to  >  0 
p{wjHi)  =  a;)^*"*  ; 

■  Ofott'i  Jo 

exp(-“«»^/£Ko)exp[“-(in“a;)^*/^il^®  (*^7) 


The  parameters  bo  and  ao  pertain  to  the  clutter  while  6i  and  ai  pertain  to  the  target. 

The  integral  in  equation  77  cannot  be  expressed  in  closed  form,  making  it 
difficult  to  determine  the  P^.  The  integral  can  be  approximated  numerically  using 
techniques  such  as  the  trapezoidal  rule  or  Simpson’s  rule  (11:368-390).  The  pdf  of 
equation  77  can  then  be  evaluated  for  a  range  of  values  of  w  and  the  distribution 
corresponding  to  the  pdf  of  77  becomes  a  discrete  distribution. 


Proceeding  using  a  binary  decision  approach,  the  likelihood  ratio,  A(w),  can 
be  written  as 

p{wlHt) 


A{w)  — 


p{w/Ho) 


(78) 


and  the  likelihood  ratio  test  becomes 


A(ti;)  = 


where  A  is  determined  by  the  particular  decision  criteria  that  is  employed.  Using  the 
Neyman'Pearson  decision  criteria  which  requires  that  A  be  determined  according  to 
a  preselected  value  of  Pja.  the  detection  threshold,  u/y,  can  be  found  by  equation  11 
and  A  becomes  (13:28-38) 

A  =  A(u;t)  (8C) 

When  dealing  with  discrete  distributions,  the  Pd  can  be  determined  by  sum¬ 
ming  up  p(wi/Hi)  for  all  i  sucl*  that  A(iw,)  >  A(iU3-)  (17).  That  is 

Pd^  ^p{v>i/IJi)  such  that  A(ufi}  >  A(w3<)  (^1) 

i 

where  i)(wi/Ni )  is  determined  from  the  numerical  evaluation  of  equation  77,  p(wi/IIo) 
is  determined  from  equation  76,  A(u;j)  is  deterimned  from  equation  VS,  and  Wi  can 


be  any  value  between  0  and  infinity.  Equation  81  can  also  be  written  as 

Pd  =  l  —  suck  that  A{wi)  <  A{wt)  (82) 

t 

The  details  of  performing  the  numerical  integration,  developing  the  discrete  distri¬ 
bution,  finding  the  relationship  between  P/a,  and  SCB,,  and  applying  the  results 
to  the  radar  range  equation  are  left  to  future  researchers. 

5.8.2  Other  Recomendations.  Future  research  should  examine  the  target  de¬ 
tection  problem  with  more  than  one  observation  of  the  returned  pulse.  That  is, 
instead  of  using  single-hit  detection,  utilize  pulse  integration  in  the  detection  scheme 
and  examine  the  problem  with  multiple  pulse  observations. 

Also,  the  problem  could  be  extended  by  determining  the  detection  characteris¬ 
tics  when  other  types  of  radar  are  used.  A  comparison  of  the  detection  characteristics 
could  be  made  for  pulse,  moving  target  indicator  (MTI),  pulse  doppler,  and  contin¬ 
uous  wave  (CW)  radais. 

5.3  Summary 

By  using  theexpressions  for  Pd  developed  in  Chapter  3  and  the  algorithm  for 
evaluating  these  expressions  presented  in  Chapter  4,  the  detection  diaracteristics 
of  a  steady  Target  in  WeibuU  clutter  can  be  determined.  Pd  can  be  found  to  a 
desired  precision  for  varying  values  of  SCll,  Pja^  <^*id  b  and  plots  cair  be 
developed  which  show  the  relationship  between  Pj,  P/o,  SCR,  and  the  radar  detection 
range.  Finally,  the  radai'  detection  range  of  steady  targets  in  WeibuU  clutter  can  be 
predicted  by  exaiuiniug  the  plots  of  P<i  versus  H/Po> 
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*  TRis  progran  BUS  developed  in  order  to  generate  data 

*  to  be  ueed  to  plot  the  probabilty  of  detection  versus 

*  the  signal~to~cltttter  ratio  and  the  probability  oi 

*  detection  versus  the  datection  range  for  various  '  ^ 

.*  values  of  probability  of  false  ala:m  and  Weibull  \ 

*  skeuness  paraBoteys. 

■ifi 

dottblo  precinitm  3cy,pfa,bR,alpha»b«ta<,lna*lnbi^a»b»pd>di, 
.+  atw!l,,ic,r^e 

-■real  y,arg  '  -  y 
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♦  The  tfeihuU  shaiKieBS  pariiBeter  <bw)  the  probabilty  o#  * 
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;  ..print,#,  ‘ESm  Te&UEIBULL  SKEWIESS  PAiy^ 

■■  -  -read  *,  bn  ■  . .  . .  ' 

print  ♦,  *EK?KR  TUB  PBQBABILJ^Y  OP  !PALSB  AUR^*  ' 
rend:  *,  pia,  '  "  ■  , 

SOD  print  »,  ^ESTER  A  i  IF  TEE  eRAaihO  lOJglB.lS  ,LAK^^ 
priiit  *,  ♦<beh«.¥idth'-Xiiait«d  .case)*  : 

'  *,  .'*E»tB&-A  2  IP  THE  GSA«SJ8G  ,AlfSLB.;TS.  SiaLU 

print  ♦,  .Kpulee-lengfch-liaited.oaga^e^  ; 

■..-'■...read*,  angle  ■  /;  ■■■-.  '  ■  .c.yy. 

vy'-if/'Caag.!*' ...be.  i  .and,  angle- -»«•  2),.  then 

y  :  --V. 

.yv end  if..-. ■■  '■■y"^';  -•  ^ 
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if  (pfa  .ge.  ld-7)  scr  =  O.dO 
if  (pfa  .It.  Id-7)  scr  =  8.d0 
end  if 

if  (bv  .ge.  1.8d0  .and.  b«  .It.  2.d0)  then 
if  (pfa  .ge.  l.d-7)  scr  =  O.dO 

if  (pfa  .ge.  l.d-il  .and.  pfa  .It.  Id-7)  scr  =  10. dO 
if  (pfa  .It.  l.d-11)  scr  =  14. dO 
end  if 

if  (bW  .ge.  1.6d0  .and.  b9  .It.  1.8d0)  then 
if  (pfa  .ge.  l.d-5)  scr  ~  O.dO 

if  (pfa  .ge.  l.d-9  .and.  pfa  .It.  l.d-5)  scr  =  12. dO 
if  (pfa  .It.  l.d-8)  scr  =  17. dO 
end  if 

if  (bw  .ge.  1.4d0  .and.  bv  .It.  1.6d0)  then 
if  (pfa  .ge,  l.d-4)  scr  =  O.dO 

if  (pfa  .ge,.  l.d-7  .and.  pfa  .It.  l.d-4)  scr  =  14. dO 
if  (pfa  .ge.  l.d-10  .and.  pfa  .It. ' l.d-7)  scr  =  19.d0 
if  (pfa  .It.  l,d-}0)  scr  =  22. dO 
end  if 

if  (bv  .ge.  1.2d0  .and.  bv  .It.  1.4d0)  then 
if  (pfa  .ge.  l.d-3)  scr  =  O.dO 

if  (pfa  .ge.  l.d-6  .and.  pfa  .It.  l,d-3)  scr  =  17. dO 

if  (pfa  .ge.  l.d-9  .and.  pfa  .It.  l.d-6)  scr  =  23. dO 

if  (pfa  .ge.  l.d-12  .and.  pfa  .It.  l.d-9)  scr  =  26. dO 
if  (pfa  .It,  l.d-12)  scr  =  28.d0  ^ 

end  if 

if  (bw  .ge.  l.dO  .and.  bv  .It.  1.2dO)  then 
if  (pfa  .go.  i.d-2)  scr  =  O.dO 

if  (pfa  .ge.  l.d-4  .and.  pfa  .It.  J,d-2)  scr  =  20,d0 

if  (pfa  .ge.  l.d-6  .and.  pfa  .It.  l.d-4)  scr  °  25.d0 

if  (pfa  .ge.  l.d-0  .and.  pfa  .It,  l.d-8)  scr  -  29. dO 

if  (pfa  .ge.  l.d-11  .and.  pfa  .It.  l.d-0)  scr  =  32. dO 
Ji  (pfa  .ge.  l,d-13  .and.  pfa  .It.  l.d-11)  scr  =  34. dO 
if  (pfa  .It.  l.d-13)  scr  »  3S.dO 
end  if 

if  (bv  .8d0  .and.  bw  .It.  l.dO)  than 
’;if  (pfa  ge.  l.d-2)  scr  =  8.d0 

if  (pfa  .go.  l.d-4  .and.  pfa  .It.  l.d-2)  set  ^  27.40 

if  (pfa  .go.  l.d-8  .and.  pfa  .It.  l.d-4)  scr '»  33. dO 

.if  (pfa  l.d-8  .and.  pfa  ,lt.  l.d-6)  scr  »  37. dO  ^ 

if  (pfa  .go.  l.d-9  .and.  pfa  .It.  l.d-8)  scr  »  40. dO 

.  if  (pfs|,  ,ge.  l.d-lO  .and.' pfa  .It.  l.d-9)  scr  «  4i.d0 
if  (pfa  .ge.  l.d'^11  .and.  pfa  .It.  l.d-10)  scr  e  42. dO 

if  (pfa  .ge.  l.d-ia  .and.  pfa  .It.  l.d-il)  scr  js  43. dO 

if  (pfa  .ge.'l,d-13  .and.  pfa  .It.  i.d-12)  scr  »  44.d0 
if  (pfa  , It.  i.d-13)  ac-'  c4S.dO 
.end  if  . : 

if  (bv  :.ge.  .6d0  .and.  bv  .It.  .8d0)  then 
if  (pfa  .ge,  l.d-i),  acr  *  19. dO 

if  (pfa.r.ge.  l.d-2  .aitd.  pfa'.lt.  l,d-4)  scr  =  31. dO 

(pfa  .gs.  l,d-3  .and*,  pfa '.4t.  l.d-2)  scr  »  37. dO 


60 


at 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


if 

(pfa 

.ge. 

1. 

,d-4 

.and.  pfa 

.It. 

l.d-3) 

if 

(pfa 

•ge. 

1, 

.d-6 

.and.  pfa 

.It. 

l.d-4) 

if 

(pfa 

•  ge. 

1. 

,d-6 

.and.  pfa 

.It. 

l.d-6) 

if 

(pfa 

•ge. 

1. 

,d-7 

.and.  pfa 

.It. 

l.d-6) 

if 

(pfa 

.ge. 

1. 

,d-8 

.and.  pfa 

.It. 

l.d-7) 

if 

(pfa 

•ge. 

1. 

,d-9 

.and.  pfa 

.It. 

l.d-8) 

if 

(pfa 

.It. 

1. 

.d-9) 

scr  3  55, 

dO 

end  il 

if  (bv  .It.  .6d0)  then 

if  (pfa  .ge.  l.d-1)  scr  =  22. dO 
if  (pfa  .ge.  l.d-2  .and.  pfa  .It.  l.d-1) 

if  (pfa  .ge.  l.d-3  .and.  pfa  .It.  l.d-2) 

if  (pfa  .ge.  l.d-4  .and.  pfa  .It.  l.d-3) 

if  (pfa  .ge.  l.d-6  .and.  pfa  .It.  l.d-4) 

if  (pfa  .ge.  l.d-6  .and.  pfa  .It.  l.d-5) 

if  (pfa  .ge.  l.d-7  .and.  pfa  .It.  i.d-6) 

if  (pfa  .ge.  l.d-8  .and.  pfa  .It.  l,d-7) 

if  (pfa  .ge.  l.d-9  .and.  pfa  .It.  l.d-8) 

if  (pfa  .It.  l.d-9)  scr  3  66.4d0 
end  if 

*  Begin  loop  to  determine  the  probability  of 

*  for  current  value  of  scr. 

do  30  j  3  1,10000 

*  Convert  scr  from  dB  value  * 

scr  3  I0.d0*o>(8cr/10.d0) 


scr  =  41. dO 
scr  3  44. do 
scr  3  47. dO 
scr  3  50. dO 
scr  3  62. dO 
scr  3  53.5d0 


scr  3  35. dO 
scr  3  45. dO 
scr  3  50. dO 
scr  3  54. dO 
scr  3  67.6d0 
scr  3  60.2d0 
scr  3  62.5d0 
scr  3  64.6d0 


detection  (pd)  * 

e 


*  Compute  the  arguments  of  Marcum's  Q  function  (alpha  and  beta)  * 

*  lr<»  equation  19.  * 


alpha  =  (s«»')'**(l.d0/2.d0) 

beta  3  (2  *  dlog(l.d0/pfa))'*<»(l.d0/bw) 

C 

C  *  The  probability  of  detection  is  the  Q  function  evaluated  at  » 

C  .  •»  these  values  of  alpha  and  beta.  Compute  the  initial  value  * 

C  ;  *  of  pd  according  to  equations  62  through  66  e 

C 

Ina  a  -(alphaee2.d0)/2.d0 
Inb  3  -(bBta**2.dO)/2.dO 
a  0  desp(lna) 
b  3  dexp(lnb) 
pd  3  a  *  b 
C 

C  *  Begin  loop  to  determine  pd  according  to  the  recursive  « 

*  relationship  given  by  equation  61.  e 

.  do  10  i  3  1,2147483646 

C  .  '  '  ■ 
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C  di  is  a  double  precision  value  oi  tbe  integer  i  * 

c  '  :  ■ .  ‘ 

di  =  i 
C 

C  *  storl  stores  tbe  previous  value  of  pd  * 

C 

storl  =  pd 
C 

C  ’!<  pd  is  determined  according  to  equations  67  tbrougb  71  * 

C 

Ina  =  Ina  +  dlog((alplia**2.d0)/2.d0)  -  dlog(di) 

Inb  =  Inb  +  dlog((betaee2.d0}/2.d0)  -  dlog(di) 
a  =  dezp(lna) 
b  =  b  +  dexp(lnb) 
pd  =  pd  +  a  *  b 
C 

C  *  H  the  current  and  previous  values  oi  pd  satisfy  equation  72, 
C  *  then  exit  loop. 

C 

if  (pd.gt.ld-SSO.and.storl/pd.gt. .999999)  go  to  20 
10  continue 
C 

C  *  Convert  the  peah  signed-to-*noiae  ratio  (scr)  to  the  * 

C  •  signal-to-clutter  ratio  (scr)  according  to  equation  28.  * 

C  *  arg  represents  the  sirgument  of  the  gamma  function.  The  * 

C  <*>  gamma  function  is  evaluated  using  the  IHSL  library.  e 

C  e  The  result  of  the  gamma  function  is  stored  as  z.  *  ■ 

C 

20  arg  s  1,  2./bu 

X  a  gamma(arg) 

scr  °  (scr/2.d0)/(x*e(b»/2.d0)) 

C 

C  *  For  small  grazing  angles,  the  range  (range)  is  computed  * 

C  e  according  to  equation  39.  * 

C 

if  (angle  .eq.  2)  then 
range  =  l.dO/scr 
else 

*  For  large  grazing  angles,  the  range  is  computed  * 

*  according  to  equation  36.  * 

range  »  (8cr}**(*‘i.dO/2.dO} 
end  if 

*  Convert  range  and  scr  to  dB.  * 

range  a  lo  ♦  dloglO(range>  ^ 
scr  iO.dO  *  dloglO(scr) 

*  pd  will  be  the  y-axis  in  the  plots  of  the  detection  * 
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C 

C 

c 

c 


c 

c 

c 

c 


c 

c 

c 


*  characteristics.  To  obtain  plots  similar  to  those  * 

*  produced  in  radeur  literature,  the  y-axis  ueis  * 

*  adjusted  as  follows:  * 

if  (pd  .gt.  .0001  .and.  pd  .It.  .9999)  then 


if 

(pd 

.le.  .001)  y  =  (loglO(pd)  +  4.)  * 

8.64 

if 

(pd 

.gt.  .001  .and.  pd  .le.  .01) 

+ 

y  = 

(loglO(pd)  +  3.)  *  (18.87  -  8.54) 

+  8.64 

if 

(pd 

.gt.  .01  .and.  pd  .le.  .1) 

+ 

y  = 

(loglO(pd)  +  2.)  ♦  (32.82  -  18.87)  +  18.87 

if 

(pd 

.gt.  .1  .and.  pd  .le.  .2) 

+ 

y  = 

(Iogl0(pd)+l.)*(l./logl0(2.))*(38 

.87-32. 82)+32. 82 

if 

(pd 

.gt.  .2  .and.  pd  .le.  .3) 

+ 

y  = 

10.  *  (pd  -  .2)  *  (43.23  -  38.87) 

+  38.87 

if 

(pd 

.gt.  .3  .and.  pd  .le.  .7) 

+ 

y  = 

(6./2.)  *  (pd  -  .3)  *  (66.94  -  43 

.23)  43.23 

if 

(pd 

.gt.  .7  .and.  pd  .le.  .8) 

+ 

y  = 

10.  *  (pd  -  .7)  *  (61.13  -  66.94) 

*  66.84 

if 

(pd 

.gt.  .8  .and,  pd  .le.  .9) 

+ 

y  = 

(Iogl0(l.-pd)+l.)*(l./logl0(2.))e 

(61, 13-67. 26)+67.; 

if 

(pd 

,gt.  .9  .and.  pd  .le.  .99) 

+ 

y  = 

(loglOCl.  -  pd)  +  2.)  ♦  (67.26  - 

81.37)  +81.37 

if 

(pd 

.gt.  .99  .and.  pd  .le.  .999) 

+ 

y  = 

(loglOd.  -  pd)  +  3.)  •  (81.37  - 

91.63)  +  91.53 

if 

(pd 

.gt.  .999) 

.+ 

y  = 

(loglOd.  -  pd)  +  4.)  *  (91.53  - 

100.)  +  100. 

*  If 

pd  is  between  .0001  and  .9999,  output  scr, 

range,  pd,  and  y 

*  which  can  be  used  to  develop  plots.  * 

print  '(fl0.4,fi0.4,fl0.6,fi0.4)*,  scr,  range,  pd,  y 
end  if 

*  Stop  when  pd  is  greater  than  ,9999  * 

if  (pd  .gt.  .9999)  stop 


e  Convert  the  signal-to~clutt6r  ratio  (scr)  bach  to  the 

*  peak  signol'-'to-noise  ratio  (scr).  * 

scr  =  10.d0**(8cr/10.d0) 
scr  a  scr  *  a.dO  *  x**(bw/2.d0) 
scr  a  10. dO  e  dloglO(scr} 

•  Increment  scr.  * 

scr  »  scr  +  .01 
e  hepeat  calculation^.  * 

30  continue 

(i3  ^ 


Appendix  B.  Weibull  Probability  Density  Fmction 

The  Weibull  pdf  described  by  equation  7  has  two  parameters,  a  and  b.  The 
parameter  a  is  related  to  the  median  of  the  Weibull  distribution,  Wm,  by  equation  8, 
while  b  determines  the  skewness  of  the  distribution.  Figure  31  is  a  plot  of  the  Weibull 
pdf  for  several  values  of  b  with  tOm  =  !• 
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Figure  31.  Weibull  pdf 
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